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D 

Abstract. Griiner put forward a single particle model of charge-density wave, which is a typ- 
ical nonlinear differential equation, and also a mathematical model of pendulum. This Letter 
analyzes the solution of equation by the rotated vector fields theory, providing the relation be- 
tween the applied field E and the periodic solution, and the conclusion that the critical value 
of E for the periodic solution is fixed in the over-damped situation. With these conclusions, 
it derives the formulae of nonlinear conductivity, narrow-band noise, which are consistent with 
the empirical ones given by Fleming. 
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1. Introduction 

Charge-density waves (CDW) is a special phenomenon of electronic condensate in the low 
dimensional materials and observed in the High Temperature Superconductor materials as 
well [1]. Many special features associated with CDW, for example, the nonlinear conductiv- 
ity, narrow-band noise, have been attracting great attention [2-5]. Describing the nonlinear 
conductivity and narrow-band noise, Griiner and his collaborators proposed a single particle 
model [6] , of which the equation is 

where (3 — E/Eq with Eq = (X/2Tr)(muj 2 /e) is constant and E is an electric field applied, 
r is the friction coefficient. Equation ([!]) is called Gruner's equation in the CDW theory as 
well as is the mathematical model of driven, damped pendulum. This is a typical nonlinear 
differential equation which is non-integrable. In the over-damped situation, Griiner gave an 
approximate solution with neglect of the inertial term [6]. The derived conductivity, i.e. <j(E) = 
(Jo a/ E 2 — Eq , has an excellent behavior served as the exponent law given in Fleming's paper in 
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the region of high field, but this neglect obviously loses some characteristics of original equation, 
such that it leads to a divergence of the derivative at the threshold field, a intrinsic difference 
from the experiment result [7]. 

There are two approaches in the research of nonlinear differential equation. The first one is 
a classical method which attempts to express the solutions of differential equations into primary 
functions named closing, or into power series form. This method is difficult to develop further. 
The second is so-called " the vector field theory of differential equation" , originated by Henri 
Poincare [8], which treats the solutions as integral curves in the phase space and derives the 
qualitative properties of the solution geometrically. Many scientists have studied in this theory's 
area, such as, G. F. Duff proposed the rotated vector fields in 1953 [9], subsequently, G. Scifcrt 
et al. developed it to the general rotated vector fields [10-12]. 

With general rotated vector fields, this Letter provides some conclusions on the solutions 
of equation (JT|), that is, the relation between (3 and the periodic solution, and that, with T big 
enough (over damping), the critical value (3q wm remain unchanged, i.e., (3$ = 1. Combining 
these conclusions with the multiple segments model advanced by Portis [13], it derives the same 
formulae as the empirical ones given by Fleming [7]. 

2. The periodic solution of the equation 

On the cj) — z phase plane, Equation (111 takes the form of vector field: 



where r > 0, j3 > 0, whose solutions correspond to the trajectories on the plane. Each 
trajectory, has a direction running as time goes forward, to which the tangents become vectors 
which constitute the vector field. When the vectors are rotating owing to the change of some 
parameter of equation, they constitute the rotated vector field. A point satisfying ^ = and 

= is called a "singularity" . A field, in which the singularities happen to move accompanying 
the rotation of vectors, is called the general rotated vector field. In the theory of differential 
equation, owing to the uniqueness of the solution, the trajectories do not intersect each other 
except at singularities [14]. This is very useful property in the following analysis. 

While (3 > 1, there exist no singularities. While < (3 < 1, there exist singularities on the 
$-axis with coordinates denoted by (<p n ,0) where 



which are categorized into two classes, one denoted by Ak(4>2k-i, 0), the other denoted by 
-Sfc(02fc,O), k £ Z. While T > 0, Bk(4>2k,0) is called a "focus", which is stable, nearby which 
all the trajectories spiral into it as time goes forward, as shown in Figure [l](a). While T = 0, 
it is called a "center", which is not the subject of the current analysis. Ak((f>2k-i, 0) is called 
a "saddle", nearby which there are four trajectories approaching it along two separatrices, as 
shown in Figure [11(b). These singularities are key to analyze the equation (pj. 




(2) 




(3) 
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Figure 1: The plot of singularities: (a) focus, (b) saddle 



For convenience, let W((f>, z) = z, Q(4>, z) = [3 — sine/) — Tz. The angle between the vector 
and the $-axis is denoted by 6, where 

Clearly, the vector field constructs a rotated vector field when parameter T or /3 changes. It 
is different that the singularities do not move as parameter T changes but move as parameter 
(3 does, so it is a general rotated vector field with respect to (3. The distribution of vectors is 
determined by W((j),z) and Q(4>,z), both with a period of 2n along the $-axis, so the vectors 
distribute with the same period and the plane can be rolled up into a cylindrical surface, that is 
called a cylindrical system in mathematics. Thus it is enough for analysis that we just discuss 
the interval [— tt— </>o, 7r— <Pq] without special declaiming, whereas the others can be extended with 
periods along the <I>-axis. In the cylindrical system, the trajectory, which connects successive 
periodic points, will become a cycle on the cylindrical surface, as shown in Figure [2] which is a 
periodic solution. It is found that equation ([2| just possesses such a periodic solution in some 
condition. 




Figure 2: Cycle on the cylindrical surface 

It is very difficult to tackle the equation ([2| directly, but in the case where T = and (3 = 0, 

it offers a simple equation, i.e., — — +sin0 = 0, which is called sine-Gordon's equation, for 

at 2 
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which it is easy to give a solution, 



z = ^2(cos^r+I). (5) 
Its corresponding trajectory Lq connects saddle point Aq with saddle point A\, as shown in 
Figure [3j Point A\ is a periodic point of Aq, so trajectory L is a cycle on the cylindrical 
surface, a periodic solution. As connecting two successive saddles, it is called a critical periodic 
solution in mathematics. 




Figure 3: Trajectories as (3 = 



For which, let T increase from zero to a positive, by the analysis of rotated vector field, it is 
proved that trajectory L will roll down and become trajectory R as shown in Figure [3] While 
r > 0, let P increase from zero to a positive value, it is proved that trajectory R will stretch 
up and singularity A\ moves toward left as well, such that trajectory R enters singularity A\ 
again at a certain positive value of (3 which is no greater than 1 [15]. Since point A\ is the 
successive periodic point of Aq, trajectory R connecting these two saddles will become a cycle 
on the cylindrical surface, a critical periodic solution, as shown in Figure [2] 

In other words, equation ^ must possess a value (3q for arbitrary T > 0, while (3 = (3q, 
equation ^ has a critical periodic solution which connects saddle A with saddle A\. By the 
analysis of rotated vector field, it can also be proved that, while (3 > 0o, equation ^ possesses 
a unique and positive periodic solution, which is stable as well; while f3 < (3q, equation ^ 
possesses no periodic solution. Through analysis, it shows that (3q is unique for a given T > 0. 

According to equation it is easy to obtain 

dz /3-sin^ 

T4> = ——- T - (6) 

On the base of equality , it is easy to obtain Theorem 1 as follows. 

Theorem 1. If a periodic solution z = z{4>) exists for the equation Q), i.e., z(tf) = z(<j)+ 2tt), 
4> S (— oo,+oo), then it must satisfy 

z(0)d0=-^. (7) 
4> L 
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The proof of Theorem 1 is very simple. According to equality (|6j, integrating from <fr to 
(j> + 2n and noticing the periodicity of z(cj>), it is easy to obtained equality 

According to Theorem 1, it is easy to prove that the periodic solution is unique and positive. 
In the proof, it just need to notice the property that the trajectories do not intersect each others. 

According to the following reason, it is easy to prove that the periodic solution is stable. 
The periodic solution is a cycle on the cylindrical surface, of which the characteristic exponent 
is — r < 0, therefore it is a stable limit-cycle [14], a stable periodic solution. 

Supposing (3o(Ti) = 1 and T2 > Ti, it can be proved that, Po(^2) = 1, that is to say, the 
critical value /3 equals 1 all the time if T is bigger enough. Interestingly, the slope of -R(r 2 ) at 
point A\ is always 0, as shown in Figure [4] 




3^ -4 

2 



-2 -1 




Figure 4: The critical periodic trajectories as T-2 > Ti 



3. Interpretation of the periodic solution in physics 

According to the conclusion above, it follows that while (3 > /3q, Equation ([2| must possess 
a unique, stable and positive periodic solution, where a stable trajectory means that all those 
nearby will approach it as time goes forward. It indicates that, whatever the initial state it is, 
the final state will be in a periodic motion toward the positive direction, as shown in Figure [5] 
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Figure 5: A stable, positive periodic solution 
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While (3 < Po, Equation ([2| possesses no periodic solution. In fact, all trajectories on the 
phase plane will approach the saddles or foci. If there are some disturbances, they will deviate 
from the saddles, which are unstable; finally, all the trajectories will end at the foci, which are 
stable, as shown in Figure [6] In the physics, whatever the initial state it is, the system will stay 
at equilibrium points (foci) in the stable state. 
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Figure 6: Nonperiodic trajectories 

Supposing a minimum value r m i n such that 0o(T m m) = 1, as above it follows that, while 
r > r min , the critical value (3 will be 1 all the time. M. Urabe has given the definite value of 
it by numerical calculations [16], that is, r m ; n w 1.193. According to the experimental values 
given in the article [13], uj « 2tt x 210(MHz), t = 1.3 x 10~ 10 (sec), T = (w t)~\ it derives 
r = 5.83. Clearly, in the practical case, T is greater than r m j n so far so that /3q = 1. Because 
[3 = E/Eq, it follows that, while E > Eq, Equation (JTJ possesses a unique periodic solution 
which satisfies Theorem 1, i.e., 

2nE r2 * 

Letf =(f) c +(f) o , where 



r^o J ^ (8) 



dtj c TE ' 19 > 

(^f) ^ s um f° rrn velocity of moving as an ensemble and contributing to the direct current 

conductivity, (^fj IS the modulated part producing the narrow-band noise. As the differ- 
ence of velocity between segments results the pressing and stretching each other, the constant 
component in (^fj > i- e - 

1 



dt;/-r< < u > 

/ a I 

is turned into energy damage so that no contribution to the current. 
4. Application of the multiple segments model 

According to the opinions given by Portis [13], the CDWs are pinned by both weak and 
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strong impurities. The strong impurities, separating the CDW into small segments, are ran- 
domly distributed leading to a distribution of separations lc 

P(y = ~e-^>. (12) 

where (Zj) is the mean separation between strong impurities. Each segment is pinned by 
one strong impurity such that the pinning potential per length acting by strong impurities 
is ((Zj) /Zj) sin </>. The weak impurities are present in higher concentration such that treated as 
a constant per length on the average and equivalent to a negative field Ex- Thus, the equation 
for the Zith segment of CDW is expected to be of the form: 

where the friction coefficient T arises from thermodynamic damping. On the analogy of Equa- 
tion ([I]), it derives that when the field E is large than a critical field Ei, this segment will slip, 
where 

E i = E T + ^-E . (14) 

w 

According to equality (14 1, it derives that, at a given electric field E there is a critical CDW 

length Ie such that CDWs longer than Ie are slipping and CDWs shorter than Ie stay at the 

equilibrium points (i.e. foci) where 

Ie = (U) t^V- (15) 
As all CDWs connect together, integrating Equation (13), it obtains 

+r^}W(W^ = [ + °°{^(E-E T ) - (k) S m<j>}P(kA. (16) 







di* dt> ^"(li) J lE L iV J/ w/ " J W/ (Z, 



Completing the integral of Equation (16 1, it obtains 

^f + T^ + e e - e t sin0= — [E + (E-E T )]e e ~ e t . (17) 
Let vcdw denotes the velocity by which the ensemble of CDWs move together, it follows that, 

Xu ( d<j)\ 

vcdw = -a-[-u) ■ ( 18 ) 



2?r \dt, o 

The factors Xujo(2tt)^ 1 appear above because Equation (17 1 is dimensionless and dcf>/dt is a 
phase velocity. Comparing Equation (17 1 with (|T|), it obtains 

vcdw = ~ ■ ~(E - E T )e~^ . (19) 

ZTT 1 CjQ 

Thus, the conductivity contributed by CDW is 

nevcDW H E T - ®l 
°dc = ^ =cr (l-— )e e - e t, (20) 

where <7o is a measure of the CDW conductivity and n is the concentration of condensed carriers 
(i.e. CDW). This formula of conductivity is the same as the part contributed by CDW in the 
Fleming's empirical ones [7]. Meanwhile, the fundamental frequency of narrow-band noise is 

, vcdw Q t (oU 
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where Q is the wave vector of CDW, Jdc is the direct current density contributed by CDW. 
The formula (21 1, that the fundamental frequency is proportional to the current, is consistent 
with the result in the article [7], which is postulated on experimental ground. 



5. Discussion of Conclusions 

The driven damped pendulum is a essential model of nonlinear science [17, 18], which de- 
velops many manifestations including chaos and appears in many physical subjects such as the 
Josephson junctions [19-21]. Delicately, the CDW and the Josephson junction, revealing the 
chaotic behaviors of conducting [22,23], are both well described by the driven damped pendulum 
model [24] . By a rigorous mathematical analysis are derived these conclusions of the pendulum 
equation applied in [15], which provide us a more direct relation between the pendulum model 
and the behavior of the CDW. Such a simple model of pendulums indicates us an interesting 
relations of these physical subjects. 
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